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Abstract 

We study the abelian sandpile model on decorated one dimensional chains. 
We determine the structure and the asymptotic form of distribution of 
avalanche-sizes in these models, and show that these differ qualitatively from 
the behavior on a simple linear chain. We find that the probability distri- 
bution of the total number of topplings s on a finite system of size L is not 
described by a simple finite size scaling form, but by a linear combination of 
two simple scaling forms ProbL{s) = xf^iz) + T^f'^^J?)-' -^^'^ large L, where 
fi and /2 are some scaling functions of one argument. 
PACS numbers: 05.40.-hj, 05.70.Jk, 05.70.Ln 
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In recent years there has been a lot of interest in the systems showing self-organized 
criticality (SOC) However the precise conditions under which the steady state of 

a driven system shows critical (long range) correlations are not well understood for non- 
conservative systems [^|-|^. In the case of systems with conservation laws ^j^, for example 
in sandpile models with local conservation of sand, it is easily shown that the average number 
of topplings in an avalanche diverges as a power of the system size [^0- This, however, is 
not sufficient to ensure criticality, if criticality is defined as the existence of power law tails 
in the distribution of avalanche sizes fll] . 



Lacking a general theory, most studies of SOC depend upon numerical simulations for 
evidence of criticality. To incorporate the effect of finite size cut-offs, it is usual to fit data 
to a finite-size scaling form in which the critical exponents of the infinite system appear 
as parameters. However, on the basis of extensive numerical studies of one dimensional 
sandpile automata, Kadanoff and coworkers |[T^ , p!3[| have argued that there is more than one 



characteristic length scale in many of these models. Consequently, a simple finite size scaling 
does not describe the statistics of avalanches, and a more general 'multifractal' scaling form 
seems necessary. 

As the finite-size scaling assumption based on a single scaling form is widely used in the 
studies of SOC, it seems desirable to test it in a simple analytically tractable model. This 
we do in this Rapid Communication for a class of one dimensional Abelian Sandpile Models 
(ASM). We find that for large L, the distribution function of duration t of an avalanche, and 
the number of distinct sites toppled Sd in our model do have a simple scaling form. However, 
the distribution function of total number of topplings s, and of the maximum number of 
topplings Tie at any one site does not have a simple scaling form, but a more complicated 
linear combination of two simple scaling forms (LC2SSF) 

ProbLiX) = L-^^hiXL--^) + L-f'hiXL-"') , (1) 

for large L, where Pi = Ui = and /?2 = = 1 for X = ric and /?i = z/i = 1 and = z/2 = 2 
for X = s, and /i and /2 are scaling functions, different for X = s and X = ric. We also 



find that this behaviour is quite robust and does not depend on the choice of the unit cell, 
but in general the function fi and /2 are not universal. 



The ASM on a simple linear chain has been studied earlier by Bak et al |T^, and in 



more detail by Ruelle and Sen We consider ASM on one dimensional chains formed 
by joining a single type of unit cells (see Fig. 1). Such decorated chains are the simplest 
generalization of the linear chain. We have studied two cases in detail. Case A is a chain 
of alternating double and single bonds. Case B is a chain of diamonds joined together by 
single bonds. We solve these models exactly in the large L limit, and find that the avalanche 
distribution function shows a nontrivial behaviour, different from that of the simple linear 
chain (case C). In fact the behavior of the ASM in case C is not typical of one-dimensional 
ASM's. 

The model is defined as follows: A site on the chain is denoted by a pair of indices 
where i = 1 to L labels the unit cell and j numbers a site within the unit cell. In case A j 
ranges from 1 to 2, and in case B, from 1 to 4. At each site there is an integer variable 
hij, called height of the sandpile at that site. A particle is added at a randomly selected site. 
If the height h^j is greater than a preassigned threshold height /if^ at that site it topples, 
and loses one particle to each of its neighbours. We choose /if^ to be independent of i and 
equal to the coordination number of site of type j. A toppling at a boundary site causes a 
loss of one particle from the system. The process of toppling continues until there are no 
unstable sites. After the system is stabilized a new particle is added. 

The critical steady state is easy to characterize using the general theory of ASM's 0. 
In the steady state all recurrent configurations occur with equal probability. The set of 
recurrent configurations is characterized by the burning algorithm (see 0, also |16[)- Iii 
the burning algorithm, the sites can be burnt in any order. We choose the convention that 
the burning starts from the left boundary and continues rightward as long as possible. The 
unit cell where the rightward burning stops will be called the break point. Afterwards, the 
burning is allowed to proceed leftwards from the right boundary. It is easy to see that in a 
recurrent configuration of model A, the allowed values of {ha, hi2), for i on the left of the 



break point, are (3, 3) and (3, 2). For i on the right of the break point these are (3, 3) and 
(2, 3) and at the break point these are (2, 3), (3, 1) and (1, 3). Since each doublet other than 
the break point has only two possible configurations, the entropy per site, defined as the 
logarithm of the total number of recurrent configurations divided by the number of sites, is 
finite and equal ln(2)/2 in the large L limit. For the simple linear chain, the entropy per 
site in the SOC state is zero. This fact is responsible for its non-generic behavior. 

To the left (right) of the break point the left (right) site of a doublet always has height 
3, and the probabihty of right (left) site of a doublet having height 2 and 3 is | each. The 
break point can occur at any of the L doublets with equal probability. Averaging over the 
position of the break point, this imphes that the probabilities of the left site of i th doublet 
having height 2 and 3 are i/{2L) and 1 — i/{2L) respectively. Similarly the probabilities of 
the right site of a doublet having height 2 and 3 are |(1 — i/L) and |(1 + i/L) respectively. 
Thus the average height profile in the SOC state varies linearly with i in case A, and the 
SOC state is not translationally invariant even far away from the boundaries. This feature 
is not present in case C. 

Now we describe the spread of the avalanche in model A, which again differs qualitatively 
from case C. Without loss of generality we may assume that the point where the particle 
is added to be called the source site, is to the left of the break point. Then clearly, if the 
configuration of the doublet left to it is (3, 2), the avalanche does not spread to the left and 
propagates a distance of order L upto the break point on the right. Each site affected by the 
avalanche topples only once, and the total number of topplings in an avalanche is of order 
L. Such an avalanche is said to be of type I. The probability of such avalanches is |. One 
can easily check that otherwise the avalanche propagates a distance of order L on both sides 
of the source point. In such avalanches Uc is of order L, and the total number of toppling in 
an avalanche is of order L^. Such an avalanche is said to be of type II (see Fig. 2). As the 
probability that the addition of particle will cause an avalanche is |, the fractional number 
of avalanches of type I is ^. 

The probability distributions of total number of toppling s, total number of distinct 
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sites toppled Sd, duration t, and the number of times the source site topples Uc, for type I 
avalanche can be calculated easily. It is convenient to work with the scaled variable a = |; 
and (3 = j^, such that a, /5 G [0, 1], where i and j are the position of source point and break 
point on the chain respectively. It can be easily verified that for type I avalanche s, Sd and 
t = 2{j3 — q)L. Thus the probability distribution of s/L, Sd/L and t/L for given a and /? 
for type I has a delta function at 2{j3 — a). Using the fact that a and (3 are independent 
random variables uniformly distributed between and 1, averaging over a and (3 we find for 
type I avalanches 

Pro6L(X|type I) = [1 - X/{2L)]/L , for X < 2L , (2) 

where X = s, s^, t. In type I avalanches any site topples at most once so nc = 1. 

The type II avalanches show a much complicated and interesting structure. The 
avalanche fronts, i.e. the left and right boundary sites of the active region at any time, 
do not move uniformly in time, the spreading rate depends on the local height configura- 
tion. However, for distances >> 1, one can define an average velocity. The analysis of 
these avalanches become easy using the decomposition of avalanches into waves of toppling 
proposed by Ivashkevich et al [|l^ . In each wave of toppling the source site topples only once 
and all other sites topple until they are stable. Waves of toppling propagate in exactly the 
same way as the burning front in burning algorithm. Thus a unit cell which cannot be fully 
burnt from the left (right) side stops a wave propagating towards left (right), and is modified 
so that next wave may cross it. We refer to such configurations as left (right) stoppers. The 
stoppers slow down the spreading of avalanches. Obviously the first wave propagates upto 
the break point with a velocity 1 site per time-step. To calculate the velocity towards left 
from the source point we note that (a) doublet of type (3, 3) is crossed in 2 time steps and 
(b) doublet of type (3,2) is crossed in 4 time steps because it stops the first wave approaching 
to it and it is crossed in 2 time steps by the next wave which follows after 2 time steps of the 
previous wave. Thus the average time taken by the avalanche front to cross a doublet is 3, 
which implies the average velocity is | sites per time-step. Similarly one can show that if the 
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avalanche crosses the break point on the right it will advance with an average velocity which 
is also I sites per time steps. The velocity with which avalanche front recedes backwards 
after it has hit the boundary is |. Details will be presented elsewhere. 

Since the avalanche front moves with an average velocity, it forms a polygon in the space 
time history of the avalanche (see Fig. 2) |T^. The number of sides in the polygon depend 
on the position of the source point a and break point /3 and on whether the break point is 
crossed by the avalanche or not. For example if /? > a > 5/5/6 and the break point is not 
crossed then the polygon has only four edges. If 1 — 6a > /9 > a, and the break point is 
crossed then the polygon has 6 edges. There are seven possible cases of polygons which need 
to be analysed separately. Quantities like s^, t and tlc which are proportional to the linear 
size of the polygon scale as L, and s which goes as area of the polygon scale as l? . The 
expressions of scaled variables s/L^, s^/L, tjh and can be easily evaluated in terms 
of a and /3 for each case. The probability distribution functions for given a and /3 is a sum 
of two delta functions corresponding to the cases whether the break point is crossed by the 
avalanche or not. Averaging over a and /5 we find 

Pro6i (g|type II) = ^ f' f' da d(3 Q6 {q - q,{a, (3)) (3) 

where q is the generic notation for s/L"^, Sd/L, t/L and ric/ L, Ci = 1/3 is the probability 
that the a type II avalanche crosses the break point and C2 = 2/3 is the probability that 
it does not cross the break point, and qi and q2 denote expressions of q in terms of a and 
j3 in the two cases. The full explicit expression is quite complicated and will be presented 
elsewhere. 

However, some of the important features of the distribution function can be understood 
by simple arguments. Since is the extension of the polygon along horizontal axis, s^/L is 
a linear function of a and (3 in each of the seven cases. Hence the probability distribution 
of Sd/L is a piece wise linear function. The same argument works for t and tIc also. The 
total number of topplings s is proportional to the area of the polygon. Therefore, s/L"^ is a 
quadratic function of a and (3. The probability distribution in this case is quite complicated 
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and diverges as for small {s/L'^). 

Summing over the contribution coming from type I avalanches (equation (0)) and type 
II avalanches (equation (^) we obtain the full probability distributions. Since ric and s 
scale differently for type I and type II avalanches the distributions of these quantities have 
the form given in equation (||). Other quantities like Sd and t scale as L for both types of 
avalanches. Therefore, the distribution of Sd and t have a simple scaling form. 

The treatment is easily extended to other types of unit cells also. For example in case 
B the unit cell is a diamond. In this case also, an avalanche always spreads upto the break 
point. The spread of avalanches to the other side will be either of order L or of order 1. 
Thus again, there are two types of avalanches. A detailed calculation shows that these 
occur with relative frequencies 5 : 8 on the average. While the velocities of avalanche front 
are different in this case, the probability distribution functions for both type I and type 
II avalanches have the same qualitative features irrespective of the velocities. For type I 
avalanches, t ~ ~ (/5 — a)L to order L. Thus probability distribution of Sd and t have 
same linear form as in model A, while the slope depends on the velocities. The variable 
has the probability distribution Prob{nc) ~ 2""^=. As s ~ nc{(3 — a)L, this implies that the 
scaling function /i in Eq. (|l|), is a piecewise linear function with many segments. For type 
II avalanches the space time history of active sites forms a polygon exactly like in model 
A, except that the slope of edges of the polygon depend on the velocities. Therefore the 
probability distributions have same qualitative behaviour as in model A. However the exact 
form of functions /i and /2 are not same in case A and B, and these functions are not 
universal. In case C, there are no avalanches of type I, and the simple scaling ansatz works 

In the multifractal approach one defines the function f{a) by the relation that an 
avalanche of size X = occurs with a probability which scales as L^^°'\ for large L. 
The exponent /(a) defined as lim^^oo log(-Pro6L(X)/ log(L)) is a continuous function of 
the a. For our abelian model it is easy to see from Eq. ([l|) that /(a) is a monotonically 
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decreasing piecewise linear function for X = s (see Fig. 3). We have also shown results of 
a computer simulation of the model for L = 100 for 2 x 10^ avalanches. Also shown is the 
theoretical curve using the Eq. (|1]) for L = 100 (dotted line) and L = oo (solid line). Clearly 
there is a very good agreement with simulation data. We note that The / versus a curve is 
quite similar to that obtained in and that approach to L — > oo limit is quite slow. 

As the LC2SSF involves only a finite number of unknown parameters, its use when simple 
scaling fails is preferable over the more general multifractal form. We also note that we find 
the breakdown of simple scaling without appearance of two different length scales in our 
model. 

Similar behaviour may be expected in other one dimensional models. For example, for 
ASM on L X M cylinder L>>M>>l,we expect three types of avalanches: type I and II, 
and finite avalanches of size < M, which do not ring the cylinder, and are two dimensional 
in character. This shows that a LC3SSF would describe this situation. It remains to be 
seen whether this behavior survives in higher dimensions or it is specific to one dimensional 
models. 

To summarize, we have determined an exact asymptotic finite size scaling behaviour of 
the distribution of avalanche sizes in the abelian sandpile model on a class of decorated 
one dimensional chains. We find that in these models the SOC state is not translationally 
invariant, and the probability distribution of s and ric unlike the simple linear chain is 
described by a linear combination of two simple scaling forms, and not by simple scaling 
form. 

We thank M. Barma and G. Menon for a critical reading of the manuscript. 
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Captions 

Figure 1: The one dimensional chains formed by joining (A) doublets, (B) diamonds, 
(C) single sites. 

Figure 2: The evolution of (a) type I avalanche, (b) type II avalanche in model A. The 
filled rectangle denotes a toppling event. 

Figure 3: The log — log plot of Prob{s) vs s. The solid line shows the exact asymptotic 
behaviour for L — > oo, and the dotted line shows the theoretical curve for L — 100. 
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